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We study structure of the gapless points (diabolical points) at zero magnetic field {H^ = 0) 
of single-spin models with spin anisotropics. Nontrivial appearance of diabolical points 
at finite transverse field has been studied from the view point of interference of the 
Berry phase, and related phenomena have been experimentally found in the single molecular 
magnet Fes- We study effects of the orthorhombic single- ion anisotropy E{S\ -{- S"^) and 
find a symmetry associated with the degeneracy, which provides a clear picture of the global 
structure of energy level diagram including the excited states. Moreover, we study effects 
of the higher order anisotropy C{S\_ + and find that, in contrast to the semiclassical 
limit {S oo), location of a pair annihilation of the diabolical point does not coincides 
with a point at which a pair of diabolical points appears in nonzero Hy space(bifurcation 
points). Distance between the annihilation and bifurcation points vanishes when S —> oo, 
which restores the semiclassical result. Wc obtain a complete structure of the diabolical 
points in the (C, Hx) plane. 

PACS numbers: 75.10.Jm,75.45.-Hj,75.75.-Ha,75.30.Gw,75.40.Mg,75.50.Xx 

I. INTRODUCTION 

Single molecular magnets, e.g., Mni2,Fe8 and V15, are interesting objects from both of the- 
oretical and experimental points of view in physics and chemistr}^i2ii^i^i^. Because those molecules 
consist of small number of magnetic atoms, the energy levels are discrete. There, we observe char- 
acteristics due to quantum mechanical motion of the wave function. In particular, in the high spin 
molecular magnets with an easy-axis anisotropy, such as Mni2 and Feg, a step-like magnetization 
process where suddenly changes has been observed in a sweep of the magnetic field. This 
phenomenon is understood to be attribute to the quantum tunneling between two values of M^, 
and is called resonant tunneling.— i^i^i^ii^iii. The energy level diagram as a function of the mag- 
netic field Hz consists of linear lines denoting the Zeeman energy (diabatic state) . At the crossing 
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point of the energy levels, however, they form an avoided level-crossing structure due to some 
quantum mixing interactions which cause nonzero matrix element between the crossing states. 
When the field crosses these points, the state undergoes adiabatic and nonadiabatic transitions. 
This quantum mechanical aspect of magnetization process has been studied from the view point 
of Landau-Zener-Stueckelberg mechanisn>i^ii^'i^ii^. There, the energy gap and sweeping velocity 
determine properties of the transition. By making use of this formula, determinations of the en- 
ergy gaps have been performed^^*^*^. Besides the high spin molecules, there have been also found 
various types of magnetization processes which reflect the quantum mechanical aspects of specific 
energy level diagram of the systems.— i^'^'^'^. These systems have attracted attentions also from 
view points of possible applications, for example, a basic component of a quantum computer—. 

The energy gap is understood as a tunnel splitting of the energy levels. That is, by tunneling 
between classically degenerate minima of a potential, the degeneracy is broken. The idea of quan- 
tum tunneling of magnetization was proposed by Bean and Livington^^, and the first theoretical 
description was given by ChudnovskySS. This tunneling phenomenon can be characterized by the 
instanton solution in the semiclassical treatmenta^Ii^i^. Thus, usually the ground state in finite 
quantum systems is unique. 

However, in some situation, a degeneracy can exist as has been predicted by Bogachek and 
Krive^. The point at which the energy gap vanishes is called a "diabolical point"—. It was pointed 
out that an interference of Berry phase^^ plays an important role in small magnetic particles^i^i^. 
Garg studied this phenomenon by studying destructive interference of the Berry phase by using 
the spin coherent state path integral formulation. He showed that the tunnel splitting at = 
is quenched in a single spin system of a large spin S with biaxial anisotropy of the terms {—DSl + 
E{S\ + S"^)) under nonzero transverse fields ff ^ | 36i37 gyen when Kramers' theorem is inapplicable. 
There, the tunnel splitting is found to oscillate as a function of the transverse field. That is, energy 
gaps vanish at some values of the transverse field H^- Villain and Fort studied a case of large 
spin in a weak external field limit^. They rederived Garg's result, and extended the study in the 
{Hx,Hz) plane. Kegecioglu and Garg obtained exact locations of diabolical points algebraically in 
a model Hamiltoniar>2E. 

Werensdorfer and Sessoli experimentally observed the oscillating behavior of tunnel splitting in 
the molecular magnet [Fe802(OH)i2(tacn)6]^"^ (called Feg)^. This spin system consists of eight 
Fe atoms each of which has S = 5/2 conforming a ferrimagnetic structure. The ground state 
of this molecule has the total spin S = 10^. This material is well described by a single large 
spin model. They measured tunnel splitting of this material using the Landau-Zener-Stueckelberg 
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theory. There, it is found that the number of diaboHcal points is smaller than that expected from 
S, which is called "the missing paradox" . 

Effects of the higher order anisotropy C{S^ + 5i) are also studied. Kegecioglu and Garg 
explained the missing paradox as an effect of the higher order anisotrop y^^'^^ . Bruno pointed out 
a pair annihilation of diabolical points in the (C, Hx) plane and they move to the nonzero Hy 
spaced. They discussed the case with the large S limit using spin coherent state path integral 
formulation. 

In the present paper, we point out that the mechanism of degeneracy at finite values of Hx 
can be understood from a view of a kind of parity effect in the eigenvalues of Sx which is directly 
obtained from the symmetry of the Hamiltonian of the system. This symmetry argument provides 
a clear picture of the global structure of energy level diagram including the excited states. 

Moreover, we study effects of the higher order anisotropy C{S^ + 5i) on positions of diabolical 
points in the {C,Hx) plane, and determine a complete structure of diabolical points in the plane. 
There, we find three types of pair annihilation of the diabolical points, and also find out to where 
the diabolical points move from the plane. It should be noted that, in the case of finite S, the pair 
annihilation point at finite Hx does not coincide with the point where a pair of diabolical points 
appears in nonzero Hy space (bifurcation point) in contrast to the semiclassical case {S oo)^. 
We find that the distance between the annihilation and bifurcation points vanishes when 5 — > oo. 
Namely, the semiclassical result is restored in this limit. We also study a difference in the structure 
of diabolical points for odd and even values of S, which should be related to the parity effect 
pointed in the literature^. 

This paper is organized as follows. In Sec. 2, we introduce a single spin model of single molecular 
magnets. In Sec. 3, we study symmetry of the Hamiltonian of the single-spin model in relation to 
the nontrivial degeneracy. In Sec. 4, we discuss effects of the higher order anisotropy. Finally, in 
Sec. 5, we summarize the present results. 

II. MODEL 

In this paper, we study structures of energy level diagram of a large spin model described 

by 

n = -DSl + E{Sl + si) +C{SX + Si)-H ■ S, (1) 
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where S is a spin operator with three component {Sx, Sy, Sz), H is an external magnetic field 
{Hx,Hy, Hz). The terms oi D,E and C represent the single-ion anisotropies. When D and E{< D) 



This large spin model is used to study properties of single molecular magnets such as Mni2 and 
Feg. For these molecules, the total spin S of the ground state can be regarded to be S = lO^i^. 

In particular, we study effects of the system parameters on the energy levels, and discuss the 
behavior of the diabolical points, at which the ground state is degenerate at Hz = as has 
mentioned in Introduction. Throughout the paper, we take D as a unit of energy {D = 1). 

III. SYMMETRY OF THE MODEL WITH BIAXIAL ANISOTROPY UNDER AN 

EXTERNAL FIELD 



As mentioned in Introduction, the problem of the diabolical point has been studied exten- 
sively for the model ([T|). There, the ground state degeneracy at Hz = is studied as a function of 
Hx, and found that the energy gap disappears at certain values of Hx- Generally, disappearance 
of the gap is associated with existence of a kind of symmetry. So far, the symmetry of the model 
has been discussed in the path-integral formulation, where the gap disappearance is attributed to 
a destructive interference of the Berry phase. 

In this section, we study the symmetry of the model ([T]) with C = and the magnetic field 
along X-axis: 



from a view point of explicit form of the Hamiltonian consisting of spin operators. 

Because we consider the case that the principal anisotropy axis is along the z-axis, naively we 
consider that the existence of Hx destroys the symmetry of the Hamiltonian. However, it should 
be noted that at a certain combination of D and i.e., 



are positive, the easiest axis is the z— axis {—DSl), and the hardest axis is the x— axis {2ESx)- 



A. Special symmetric point 




(2) 



E = 0.5D, 



(3) 



the Hamiltonian can be expressed as follows 



-DSl + D {Si - Si) - HxSx 
2DSl - HxSx -DS{S + l). 



(4) 
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This Hamiltonian only consists of Sx, and thus it is commutative with Sx- Therefore, this Hamil- 
tonian can be diagonahzed simultaneously with Sx, where the eigenstates are 

Sx \Mx) = Mx \Mx) , Mx = -S,-S + 1,--- ,S. (5) 

In this system, the energy levels are linear as a function of Hx, and cross each other without gap. 
Because D is positive, at Hx = the ground state is a state of Mx = 0, i.e., \Mx = 0), For S = 10, 
the ground state energy is —HOD. The first excited state is degenerate and they have Mx = ±1. 
When we increase Hx, the ground state is replaced by a state with a larger magnetization Mx + 1 
sequentially. That is, at Hx = 2, the energy level of state \Mx = 1) crosses with that of \Mx = 0), 
then \Mx = 1) becomes the ground state. Similarly, the ground state magnetization changes to 
Mx = 2, 3, • ■ ■ at Hx = 6, 10, • ■ • , respectively. In Fig. [H we depict the energy diagram of the model 
of Eq. ([2]) as a function of the field Hx ■ In Fig. [21 we plot the energy gap between the ground state 
energy (-Eg) and the first excited energy (£^1) 

AE = Ei- Eg, (6) 

by dashed lines as a function of Hx- There, we see a saw-tooth shape as shown. 




_140 I . 1 . > . > . 1 . > . J 

2.5 5 7.5 10 12.5 15 



FIG. 1: Energy diagram of the low- lying levels of the system ([T]) with S* = 10 as a function of the field H, 
for E = 0.5. 



B. General biaxial anisotropy 

Next, we consider the case with E ^ 0.5D. We set 

E = 0.bD + A. 



(7) 
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FIG. 2: Energy gap between the lowest energy and the first excited energy of the system Eq.(l2|) with = 10 
as a function of the transverse field Hx- The solid line is the case oi E = 0.485, the dashed line is the case 
of = 0.5. 



The Hamiltonian becomes 

n = no + n', (8) 

with 

n' = A {si + S^_) = 2A {Sl - S^) . (9) 

Here, the states {M^} are no more the eigenstates of the Hamiltonian Ti. The effects of the term 
Sy is expressed in terms of the raising (S^) and lowering {S~) operators of for as 

— ^ i^x'^ + '^'^''S'^, + + ^) . (10) 

This term causes the change of Mx by two. The explicit matrix element of this operator is 

{Mx = m\S^\Mx = n) 

= i (S (5 + 1) - n (n + 1))5 (5 (5 + 1) - {n + 1) (n + 2))^ 5m,n+2 
+^ {2S (S + 1) - 2n2) Sm,n 

+^ {S{S+l)-n{n- {S (5 + 1) - (n - 1) (n - 2))^ 5^,„_2. (11) 
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This term mixes the eigenstates \Mx = m) and \Mx = n) when 

\m-n\=2, (12) 

and thus it opens a gap in the crossing points with even values of |m — n\ in the energy diagram 
in Fig. [TJ In contrast, it does not open a gap between \Mx = m) and \Mx = m it 1), because 

{Mx = m\ Sl \Mx = m±l) = 0, (13) 

and, 

(M, = m\ Si \Mx = n) (M, = n| 5^ |M, = m ± 1) = (14) 

for all the possible integer values of n. 

Therefore, when the difference of the magnetization Mx between the ground state and the first 
excited state is one, the cross points in Fig. [1] remain gapless points [IS.E = 0). On the other 
hand, those of the difference two change to avoided level crossings. By this effect of S^, the 
energy diagram has a ribbon-like shape as depicted in Fig. [3l and the Hx dependence of the gap 
is smoothed as depicted in Fig. [2] by a solid curve. It should be noted that the value of E/D is 
0.082 for Feg and is much smaller for Mni2. Here we used a large value of E/D just because of 
the convenience for drawing the figure. If we use a small value of E/D, the energy difference is 
too small to see. The physical mechanism is the same irrespective of the value, and here we use a 
large value. If we decrease the value of E down to E = 0.3, the ground state and the first excited 
state almost degenerate as depicted in Fig. HI There, the energy gap IS.E has a shape which has 
often appeared in literature (Fig. [5]). 

IV. EFFECTS OF A HIGHER ORDER ANISOTROPY 

In single molecular magnets with large spins, e.g., Mni2r^'^'^ and Feg^, existence of the 
higher order anisotropic term 

n" = c{sX + st) (15) 

has been suggested. In this section, we study effects of this fourth order anisotropy. The Hamilto- 
nian without the magnetic field is 

n = -DSl + E{Sl + Sl) + n". (16) 
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FIG. 3: A ribbon like structure of energy diagram of the low-lying levels of the system Eq.© with S ~ 10 
and E = 0.485 as a function of the field Hx- 




-140 I • 1 • ' • ' • 1 • ' • J 

2.5 5 7.5 10 12.5 15 

FIG. 4: Energy diagram of the low-lying levels of the system Eq.(l2|) with = 10 and E = 0.3 as a function 
of the field Hx- The lowest energy and the first excited energy almost degenerate in this vertical axis scale. 
The second excited energy and the third excited energy are also almost degenerate. 

Here it should be noted as follows. Because S+ = Sx + iSy and S- = Sx — iSy, and 
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0.06 




FIG. 5: (a) Energy gap between the lowest energy and the first excited energy for E = 0.3 as a function 
of the field Hx- There are 10 diabolical points, (b) Energy gap between the second excited energy and the 
third excited energy for E = 0.3. In this case, there are 9 diabolical points. 



Thus, in the representation which diagonalize i.e., {|Mi,)}, it is given by 

+ St. =25^ + - [S^ + S,j.) — -S"^ [S^ + S^) — - [S^ + S^) S: 



+ Sx Sx ) + {S^ + {S^ ) Sx 

+ l{St-Sxf (18) 

which can change the value of Mx by multiples of 2. 

Therefore, nonzero components of matrix elements of the fourth term are 

{Mx = m\ n" \Mx = m) , 
{Mx = m\ n" \Mx = m±2), (19) 

and 

{Mx = m\n"\Mx = m±A) . 

Because 

{Mx = m\ n" \Mx = m ± 1) = 0, (20) 

the fact that the gap opens only at crossing points where the magnetization Mx differs by two 
maintains. 



A. Dependence on C at fixed E 

First let us study the behavior of the diabolical points on C at fixed value of E. We plot 
the change of the diabolical points in a coordinate {C,Hx) in Fig. [6l As far as \C\ is small, the 
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FIG. 6: Diabolical points between the lowest energy level and the first excited energy level on the (C, Hx) 
plane for the E ~ 0.5. The symbol (□) denotes the type I annihilation points. The symbol (o) denotes the 
type II annihilation points. 

number of diabolical points is the same as that of C = 0. However, for large \C\ cases, pairs of 
diabolical points disappear from the figure. We call this point {C,Hx) "type I an annihilation 
point" which is shown by (□) in Fig. O The pair annihilation occurs from the side of large 
when C decreases in the negative C region. In the positive side, diabolical points are drawn into 
the Hr^ axis sequentially. At the axis, the diabolical point combines with that from the negative 
Hx side, and disappears, which we call "Type II annihilation points", and denote it by (o) in Fig.[6l 
At these annihilation points, the diabolical points move to a nonzero Hy region. 

First, we show the motion of diabolical points around the type I annihilation point. In Fig. [71 
we plot the motion of diabolical points in the largest Hx values in a Hx > subspace. There, we 
find that a pair of diabolical points is created in nonzero Hy region at a point. We denote this 
point by the symbol (A). We call this point "a bifurcation point". Here, it should be noted that 
the point of the creation of the pair is not the point of the annihilation of the pair on the (C, Hx) 
plane. We find that this separation of the annihilation point and the bifurcation point exists in 
all the finite values of S. In Fig. [HI we show the case S = 2, where we find the same type of 
structure. The separation is much larger than the case of 5" = 10. 

The effect of the fourth order anisotropy has been discussed by Brunch. His argument is the 
following. There is a critical value of C = Cc where two diabolical points collide, and at this point 
the bifurcation takes place. That is, a pair of two diabolical points appears at the type I annihilation 
point. However, we find that the bifurcation point is different from the annihilation point, and 
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FIG. 7: The branch of diabohcal points between the lowest energy and first excited energy with largest 
in the case of 5 = 10 and E = 0.3. The symbol (A) denotes the bifurcation point. The symbol (□) denotes 
the annihilation point. 




FIG. 8: The branch of diabolical points between the lowest energy and first excited energy with largest 
in the case of 5 = 2 and E — 0.3. The symbol (A) denotes the bifurcation point. The symbol (□) denotes 
the annihilation point. 

appears at a larger value (smaller |C|) of C. This means that the number of diabolical points 
are not preserved on the {Hx,Hy) plane when we change C. This fact is different from Bruno's 
argument. In his arguments, the number of diabolical points on the {Hx,Hy) plane is preserved 
except at Cc On the other hand, our numerical result shows that the number of diabolical points 



12 



on the {Hx, Hy) plane can change with the value of C. Bruno's discussion is based on the large S 
limit. Thus, we study S dependence of the separation of the annihilation and bifurcation points. 

Here, we investigate structure of the diabolical points near annihilation points. In Fig. [HJ a pair 
of diabolical points near annihilation points has a parabola-like structure on the {C,Hx) plane. 
Thus, we try to fit the curve using a rotated parabola function {a'^C^-\-2abHxC+b'^H^ + cC+dHx + 
e = 0) with constants (a, b, c, d, e). The fitting is given in Fig[9l The origin of this rotated parabola 
where the diabolical point is located at (C,Hx) ~ (—0.000095,22.83468), which is indicated (O)- 
The point is not the annihilation point, and it is not the bifurcation point neither. This fact is 
indicates that the bifurcation does not occur at the origin of the parabola which is a special point 
of this figure. 

25 
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22 

-1.1-10"^ -1-10"'' -9-10"^ -8-10" 



FIG. 9: Fitting of diabolical points using a rotated parabola function. The symbols (+) denote the bare 
numerical results for the diabolical points. The symbol (A) denotes the bifurcation point. The symbol (□) 
denotes the annihilation point. The symbols (x) denote points on a rotated parabola obtained by fitting, 
and the symbol (Q) denotes the origin of the fitted parabola. (Because the scales of axes of H^. and C are 
different, the point denoted by the circle does not look like the origin.) 



Now, we study S'-dependence of the distance between the bifurcation point and the annihilation 
point. We define two quantities, 

AC = Cbif — Cann, (21) 

and 
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where Cbif and H^^^ are values of bifurcation points, and Cann and -fffnn values of annihilation 
points. 
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FIG. 10: ia.):AH^ as a function of 1/5 with E = 0.3. (b)AC as a function of 1/S with E = 0.3. 



We plot AC and AH^ as a function of 1/5", in Fig. [TOl In these figures, we find that both 
AC and AH^ rapidly decrease, when we increase S. Thus, our numerical results are consistent 
with Bruno's arguments in the large 5 limit. But, it should be noted that at finite values of S 
the bifurcation point and the annihilation point do not coincide, which indicates there exists a 
nontrivial quantum effect. 

Next, we show the motion of the diabolical points around the type II annihilation points. 
There, two diabolical points move from (C, H^, Hy = 0) to {C, Hx, Hy 7^ 0). In Fig. we show 
this motion of diabolical points in the {C, H^, Hy) space. 

As we saw above, the diabolical points disappear from the (C, Hx) plane by the pair annihilation. 
In the case that S is an odd integer, there is an odd number of diabolical points in the Hx{> 0) 
region of the {C,Hx) plane. There, the last one does not have a partner. We study how the last 
point behaves in the (C, Hx) plane. In Fig. [T2l we show behavior of diabolical points of the model 
of 5 = 3 in the (C, Hx) plane. In this case, there are three diabolical points in the region of Hx > 0. 
In Fig. [T2ya). we find the pair annihilates around C ~ —0.0039. There, the Hx value of the last 
point increases when C decreases. However, when C decreases further, it goes down and finally it 
merges to the C axis as shown in Fig. [T2rb). and merges with the partner coming from the Hx < 
region. We call this point "the type III annihilation point" . Interestingly in this case the diabolical 
points move to a nonzero H^ region (C, Hx{= 0), Hy{= 0), Hz{^ 0)). but not a nonzero Hy region 
{C^Hxi^^ 0),Hy{^ 0),Hz{= 0)) as in the other cases. We depict this behavior of diabolical points 
in Fig. [131 

In this way, all the diabolical points disappear from the {C,Hx) plane when |C| becomes large, 
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FIG. 11: Diabolical points between the lowest energy level and the first excited energy level on the 
{Hx, C, Hy) space in the case of i? = 0.3. 




-0.01 -0.0075 -0.005 -0.0025 -0.16 -0.155 -0.15 

c c 

FIG. 12: Behavior of diabolical points on the [Hx, C) plane with E = 0.3 for S* = 3 case: (a) around the 
last pair annihilates, and (b) the last one merges to the C axis {Hx = 0). 

and found three types of annihilation points. By the above studies, we figured out complete 
structure of diabolical points in the {C,Hx) plane. 



B. Dependence on E at fixed C 

So far, we studied the behavior in the (C, Hx) plane. Here let us study E dependence 
of the diabolical points. In Fig. [T^ we show diabolical points on the {Hx,E) plane for a fixed 
C(= —0.001). In Fig. [T^a). we show the case of 5 = 2, where the two diabolical points combine 
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FIG. 13: Diabolical points on {H^, C, H^) space with E = 0.3. for 5 = 3 case. 

and annihilate when E becomes small. This is a type I annihilation point. There, they move to 
nonzero Hy region. In the case of S" = 3 cases, the last one diabolical point moves to the origin 
{Hr^, E) = (0, 0) as depicted in Fig. [T^b). This is a special case of the type III annihilation point. 




1 2 3 4 01234567 

FIG. 14: Diabolical points on the {H.j,,E) plane with C = -0.001. (a) 5 = 2, and (b) 5 = 3. 

The same type behavior is found in larger spin cases (S" = 4, 5, , , ) (not shown). This observation 
indicates that the ground state for ii^ = is two fold degenerate in the odd spin cases. This is a 
degeneracy not related to Kramer's degeneracy, because S is integer. We can easily understand 
this degeneracy. For £^ = 0, a Hamiltonian is described by 

n = -DSl + C[SX + St) . (23) 

If we set C = 0, \Mz = —S) and \Mz = S) give the two fold degenerate ground state, where 
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5*2 \Mz) = Mz \Mz). For even spin cases, matrix element between the states \Mz = ±5) is nonzero 



because the difference of the magnetization (=25") is a multiple of 4, where n is an arbitrary 
integer. On the other hand, and for odd spin cases, the difference 25" is not a multiple of 4. Thus, 



Therefore, quantum tunneling between the two states does not occur, and the ground state is two 
fold degenerate in odd spin models for £' = and C 7^ cases. 



We investigated nontrivial degeneracy of eigenenergies of single molecular magnets using 
the large single spin model. In the parameter space {E,C, H^, Hy, Hz), positions of the points 
at which the eigenenergies are degenerate (diabolical points) are studied. As has been pointed 
out, the model ([T]) has diabolical points at nonzero H^. This fact seems nontrivial and has been 
studied in terms of the Berry phase in the path-integral formulation^^. We pointed out that the 
existence of diabolical points at nonzero Hx is understood from a view point of the parity effect of 
the magnetization in the x direction. 

We also studied effects of the higher order anisotropy C. For a small value of |C|, there are 
S diabolical points with positive values of H^. We studied behavior of those points when |C| 
increases. They move out from the {C,Hx) plane by pair annihilations. We found three types of 
annihilations. In the positive C case, each diabolical point moves to the C axis, and at the C axis 
it combines with the partner coming from negative region and they move to the nonzero Hy 
region. In the negative C case, the diabolical points make a pair with neighbors in the positive 
Hx region. We also found a pair creation of diabolical points in the nonzero Hy region. We should 
make emphasis that the annihilation points do not coincide with the creation (bifurcation) points 
for finite values of S. This is contrast to the case of S 00, which was studied by Bruno^. 
The asymptotic behavior in the limit S ^ 00 was studied and we found the distance between the 
annihilation and the bifurcation points decreases to zero when S increases. Thus, the argument 
of semiclassical picture is valid, but there exists an intrinsic quantum effect. In the case of odd 
integer S, one diabolical point remains unpaired and it moves to the C axis and make pair with a 
partner coming from negative Hx- In this case, we found that they move to the nonzero H^ region. 




(24) 



{Mz = S\ {Si + Siy \Mz = -S) = 0. 



(25) 



V. SUMMARY 
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